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Abstract 

We present a small x resummation for the GLAP anomalous dimension and its corre¬ 
sponding dual BFKL kernel, which includes all the available perturbative information and 
nonperturbative constraints. Specihcally, it includes all the information coming from next- 
to-leading order GLAP anomalous dimensions and BFKL kernels, from the constraints of 
momentum conservation, from renormalization-group improvement of the running coupling 
and from gluon interchange symmetry. The ensuing evolution kernel has a uniformly stable 
perturbative expansion. It is very close to the unresummed NLO GLAP kernel in most 
of the HERA kinematic region, the small x BFKL behaviour being softened by momen¬ 
tum conservation and the running of the coupling. Next-to-leading corrections are small 
thanks to the constraint of gluon interchange symmetry. This result subsumes all previous 
resummations in that it combines optimally all the information contained in them. 
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1. Introduction 


The theory of scaling violations of deep inelastic strnctnre fnnctions at small x has pro¬ 
gressed considerably over the last few years. This theory has been characterized by the 
long-standing pnzzle posed by the fact that no clear deviation of the data from the next- 
to-leading order prediction has been observed (for instance at HERA [1]), even thongh 
higher order pertnrbative corrections grow rapidly in the small x region. The recent de¬ 
termination [2] of the fnll set of three-loop (NNLO) splitting fnnctions makes a resolntion 
of this pnzzle less of an academic exercise. Indeed, resnlts obtained nsing these NNLO 
splitting fnnctions already [3] show sign of the known small x instability, for instance in 
the extraction of parton distribntions. Hence, NNLO phenomenology becomes problematic 
nnless this instability is treated in some way. 

However, the reason why snperhcially large small x corrections snm np to very mod¬ 
erate effects is now essentially nnderstood. It is the conseqnence of two effects: hrst, 
moment um conservation which provides a constraint which limits the nnbonnded growth 
of small X kernels, and second, renormalization gronp improvement which throngh the 
rnnning of the conpling softens the asymptotic small x growth. A systematic theory of 
small X resnmmation which inclndes both these effects has been developed in refs. [4-8]. In 
particnlar, we have constrncted a relatively simple analytical expression for the improved 
anomalons dimension which we argne is the best way to pnt together the information from 
LO BFKL and NLO GLAP. A shortcoming of the approach developed in these references 
is that even orders of the resnmmed pertnrbative expansion are snbject to a large model 
dependence. This is dne to the fact that the constraint of moment um conservation stabi¬ 
lizes the pertnrbative expansion of evolntion kernels to all orders, bnt only in the collinear 
region. In the “anticollinear” region, where parton virtnalities have reverse ordering (and 
thns only contribnte dne to In ^ enhancement) the instability remains. As a conseqnence, 
the stabilization of the BFKL kernel is lost in the vicinity of its minimnm, where the 
collinear and anticollinear regions weigh eqnally: the minimnm is in fact only present at 
leading order of the pertnrbative expansion. Becanse the minimnm is necessary for rnnning 
conpling resnmmation, the inclnsion of NLO BFKL resnmmation prevents it, nnless one 
makes some assnmption on the restoration of the minimnm by higher order corrections; 
this, however, introdnces a sizable model-dependence [8]. 

A similar approach to small x resnmmation which also embodies the principles of 
momentnm conservation and renormalization gronp improvement has been presented in 
refs. [9-11]. This approach differs [12] from onrs esssentially becanse it has the somewhat 
wider goal of determining the off-shell glnon Green fnnction, and it is therefore rooted in 
the BFKL eqnation. Whereas this gives for example an exact treatment of small x rnn¬ 
ning conpling corrections, it only allows a n um erical extraction of anomalons dimensions 
throngh deconvolntion from a given bonndary condition. Also, it does not natnrally lead 
to a symmetric (donble) leading log expansion of the anomalons dimension (for example, 
when Uf ^ 0 the inclnsion of snbleading GLAP terms in this approach has not yet been 
accomplished). However, in this approach it was pointed ont that the nnderlying symme¬ 
try npon the interchange of glnon virtnalities in the glnon-glnon scattering process which 
determines the BFKL kernel can be nsed for the treatment of the instability in the anti- 
collinear region. This symmetry interchanges the collinear and anticollinear regions, and 
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therefore it enables the implementation of the momentnm conservation constraint even in 
the anticollinear region. 

Im this paper, we exploit the symmetric collinear-anticollinear resnmmation [10] to 
stabilize the resnmmed and renormalization-gronp improved pertnrbative expansion of 
refs. [4-8]. This can be done by nsing the dnality relations [4,5] which connect GLAP 
and BFKL kernels in order to translate the symmetry constraints on BFKL evolntion 
into constraints on snbleading contribntions to the GLAP anomalons dimensions. Becanse 
the “donble leading” pertnrbative expansion of refs. [5] systematically inclndes leading, 
next-to-leading,.. .logs of and 4 simnltaneonsly, its symmetrization mnst be done by 
inclnsion of terms which are snbleading order by order, and this reqnires a snbtle treatment 
of donble connting terms. Also, implementation of the symmetry is nontrivial when the 
conpling rnns, becanse the choice of rnnning conpling scale may break it. Finally, after 
symmetrization, even at LO, the BFKL intercept, which determines the asymptotic small x 
behavionr, depends nonlinearly on the strong conpling, and this makes the renormalization- 
gronp improvement of refs. [7,8] rather more elaborate: it is still possible to derive an 
asymptotically exact solntion to the rnnning-conpling BFKL eqnation in terms of special 
fnnctions, bnt this now reqnires a donble expansion abont the BFKL minimnm. Once this 
is done, a stable pertnrbative expansion is obtained, and in particnlar the NLO resnmmed 
resnlt tnrns ont [13,14] to be both close to the LO resnmmed one, and also qnite close 
to the NLO hxed-order resnlt, thereby explaining the lack of experimental evidence for 
large deviations from low hxed-order predictions at small x: in fact, the resnmmed NLO 
splitting f un ction is closer to the NLO than to the NNLO hxed order one. 

In this paper we discnss the resnmmation of anomalons dimensions and splitting 
fnnctions. To make direct connection to the strnctnre fnnctions we need to also implement 
coefficient fnnctions and their resnmmation. This also entails a scheme choice: here, we 
adopt the so-called Qo scheme [15], which coincides with the standard MS scheme at hxed 
order bnt dihers from it at the resnmmed level. Issnes of scheme dependence were discnssed 
in detail in ref. [6], where resnmmed coefficient fnnctions and hts to the data were explicitly 
presented. The relation of Qo and MS after small x rnnning conpling resnmmation was 
discnssed in ref. [7]. Whereas the nse of the resnlts presented in this paper for actnal hts to 
data will reqnire the formalism presented in [6], [7], all the important small x resnmmation 
ehects are inclnded in the Qo-scheme resnlts presented here. In order to minimize issnes of 
scheme dependence, in this paper we will determine all anomalons dimensions and splitting 
fnnctions with n/ = 0. 

The paper is organized as follows. After recalling in section 2 some backgronnd resnlts 
from the resnmmation formalism of refs. [4-8], which will be needed in the seqnel, in 
section 3 we explain how the nnderlying symmetry of the BFKL kernel can be conpled to 
the donble-leading expansion of ref. [4] in order to prodnce a stable resnmmation at the 
hxed conpling level. The formalism which is reqnired to generalize the rnnning conpling 
resnmmation of ref. [7] to this symmetrized case is introdnced in section 4. Resnmmed 
anomalons dimensions and splitting fnnctions are explicitly constrncted in section 5 at 
the leading order and in section 6 at the next-to-leading order. Resnlts for anomalons 
dimensions and splitting fnnctions are presented, discnssed and compared in section 7. 
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2. Background 

The behaviour of structure functions at small x is dominated by the large eigenvalue 
of evolution in the singlet sector. Thus we consider the singlet parton density G(^, t) with 
^ = logl/x, t = \ogQ^/y?^ such that, for each moment 

G{N,t)= (2.1) 

Jo 

the associated anomalous dimension N) corresponds to the largest eigenvalue in 

the singlet sector. At large t and hxed ^ the evolution equation in A^-moment space is then 

^G(iV,t)= 7 (a,(t),iV)G(iV,t), (2.2) 

where as{t) is the running coupling. 

The perturbative anomalous dimension is known up to NNLO [2], given by 

'y{as, N) = as'yo{N) + Q;^ 7 i(iV) + Q;^ 72 (iV)- (2.3) 

The corresponding splitting function is related by a Mellin transform to 7 (as, N): 

7 ( 0 ( 3 , A^)= [ dxx^P{as,x). (2.4) 

Jo 

Small X for the splitting function corresponds to small N for the anomalous dimension: 
more precisely P ~ (l/a:)(log(l/a:))"^ corresponds to 7 ~ Even assuming that 

a leading twist description of scaling violations is still valid in some range of small x, as 
soon as x is small enough that Os^ ~ 1 , with ^ = logl/ac, all terms of order {as/x){asC)^ 
(LLx) and Os(os/a:)(os^)"^ (NLLx) which are present in the splitting functions must be 
considered in order to achieve an accuracy up to terms of order a‘^{as/x){asC)^ (NNLLx). 

These terms can in principle be derived from the knowledge of the kernel x(q(s, M) of 
the BFKL ^-evolution equation 

AG(e.Af) = x(a.,M)G(e.M). (2.5) 

which is satished at large ^ by the Mellin transform of the parton distribution 

/ OO 

dte-^^Gitt). ( 2 . 6 ) 

-CX) 


The evolution kernel in eq. (2.5) is the Mellin transform of the BFKL kernel K ^cts, 

which determines the high-s limit of the cross-section for the scattering of two gluons with 
virtualities k‘^ and into n gluons 

= (^)-"A-(„.,f^). (2.7) 
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This kernel has been compnted to NLO accnracy [16-19] 

x{ois,M) = asXo{M) + alxi{M) + .... ( 2 . 8 ) 

Notice that in glnon-glnon scattering ^ whereas in deep-inelastic scat¬ 

tering ^ ~ ln(s/(5^). The kernels that correspond to these two cases coincide at leading 
order, and beyond leading order are related in a simple way, as we shall discnss explicitly 
in sect. 3 below. 

Qnite generally, the anomalons dimension 7 eq. (2.3) and the BFKL kernel are related 
dne to the fact [20,5,6] that the solntions of the BFKL and GLAP eqnations coincide at 
leading twist if their respective evolntion kernels are related by a “dnality” relation. In 
the hxed conpling limit the dnality relations are simply given by: 

x{as,l{as,N)) = N, (2.9) 

7(«s,x(«s,M)) = M. (2.10) 

This implies that all relative corrections of order [as log l/x)"^ to the splitting f un ction can 
be derived from xo(Af), those of order as{ois^ogl/x)^ from xi(M) and so on: indeed, if 
we expand 7 ( 0 : 3 , A^) in powers of o^ at hxed ag/N 

7 ( 03 , N) = Xsi^s/N) -F 037ss(as/iV) + .. ., (2.11) 

eq. (2.9) determines 73 in terms of xo, 7 ss in terms of xo and xi and so on. 

Conversely, all relative corrections of order (oglogQ^)"^ to the BFKL kernel 

K ^ 03 , j eq. (2.7) can be derived from 7 o(A^), those of order 03(03 logQ^)"' from 71 (A^) 
and so on: expanding x(<as, M) in powers of Og at hxed ag/M 

x(o3, M) = Xs{ois/M) asXss{ois/M) -t- ..., (2.12) 

eq. ( 2 . 10 ) determines Xs in terms of 70 , Xss in terms of 70 and 71 and so on. 

As is by now well known, the early wisdom on how to implement the information from 
Xo was completely shaken by the mnch softer behavionr of the data at small x [21,22,23] 
and by the compntation of xi [16-19], which showed that the naive expansion for the 
improved anomalons dimension had a hopelessly bad behavionr, in particnlar near M = 0 
and M = 1. In refs. [6,5] we have shown that the physical reqnirement of momentnm 
conservation for anomalons dimensions hxes by dnality the valne of x{(^s, M) at M = 0 : 

X(a3,0) = l. (2.13) 

This implies that in order to cnre the instability at M = 0 it is snfhcient to ensnre that 
the standard GLAP leading log resnmmation is performed in accordance to the pertnrba- 
tive expansion of the anomalons dimension eq. (2.3), in which momentnm conservation is 
respected order by order. 

The natnral way to do this is by combining the small x resnmmation and the standard 
resnmmation of collinear singnlarities. As a resnlt one obtains a pair of ‘donble-leading’ 
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(DL) perturbative expansions of the kernel x and anomalous dimension 7 , dual to each 
other order in perturbation theory (up to subleading terms). For example, in the DL LO 
expression for x an inhnite series of terms of order (cts/M)"^ is added to the hxed-order xo 
result, and similarly at DL NLO terms of order as{c(s/M)'^. The net result is to obtain a 
regular function at M = 0 which only violates momentum by small subleading terms, and 
interpolates between the BFKL result at small x (M ~ d) and the GLAP result at large 
X (M ~ 0). 

The DL expansion kernels are 

XDL — XDL LO + OtsXDL NLO + • • • (2-14) 

with 

Tl CY 

XDLLo{as,M) = asXo{M) + Xs{^) - (2.15) 

XDLNLo{as:M) = a^Xl^M) + Xss{^) ~ ^ ~ 

The kernels Xs and Xss are derived from 7 o(lV) and 71 (A^) using the inverse duality relation 

x{as,x{o^s,M)) = M (2.17) 


and are then given by: 


7o (Xs (m)) 


M 

Cts ’ 



7i(Xs(w)) 

7o(Xs(w))’ 


(2,18) 

(2.19) 


while the numbers /o, fi and /2 are hxed so as to remove double counting between xi and 

Xss- 

The DL kernel eq. (2.14) has a stable perturbative expansion for small M, and indeed 
for all M ^ 1. However, it is still unstable near M = 1: while momentum conservation 
and duality £x x = 1 at M = 0 , near M = 1 xo behaves like 1/(1 — M) and xi behaves 
like —1/(1 — M)^. Hence, the NLO approximation for xdl has no minimum — in fact, 
the minimum is absent at any even order of the perturbative expansion eq. (2.14). This 
prevents a precise treatment of non leading terms, for which a better control of the central 
region near M = 1/2 is needed, which in turn creates a problem for the evaluation of the 
small X asymptotic behaviour, which depends on the central region. In ref. [ 6 ], however, it 
was found that phenomenology requires the all-order minimum of x to signihcantly differ 
from its low-order one, and that a successful description of the data is obtained if the 
minimum value of x is taken as a free parameter to be htted to the HERA data. 

So far we have assumed that the running of the coupling is treated perturbatively in 
the DL expansion. In particular, the running of the coupling is a leading InQ^, but sub¬ 
leading In ^ effect. As a consequence, when the coupling runs the duality relation eq. (2.9) 
is unaffected at the leading In ^ level, while at the N^LLa: it receives a perturbatively 

computable O (/locislny)^ correction. Specihcally, at NLO 


Xssiois/N) 


XliXsias/N)) XoiXsias/N))xoiXsias/N)) 

x'oiXsias/N)) 2[x'o{Xs{as/NW. 


( 2 . 20 ) 
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The running coupling correction may be viewed as an effective contribution to y: defining 


xllertiM) = Ho 


x'i(M)Xo(M) 

2Xo(M) 


( 2 . 21 ) 


7 ss eq. (2.20) is the naive (fixed coupling) dual of an effective running coupling x whose 
NLO term is supplemented by Xi°pert. (2.20). However, it is apparent that this effective 
X is singular as M —2, which implies [4] that the corresponding splitting f un ction is 
enhanced as a: —0. In fact, the contribution APs'n.{as, x) to the N"^LLx splitting function 
induced by the running-coupling duality correction to the splitting function behaves as 


APsr^{as,x) 

Ps{as,x) 


rsj 
X —>0 




( 2 . 22 ) 


The all-order resummation of these running coupling small x corrections is therefore 
needed, and it turns out [7,8] to be, on top of momentum conservation, the second crucial 
ingredient of the resummation procedure. The running coupling resummation can be 
performed [7] by taking a quadratic approximation to the BFKL kernel in the vicinity of 
its minimum. Indeed, in M space the usual running coupling as{t) becomes a differential 
operator: assuming the coupling to run in the usual way with and taking only the 
one-loop beta function into account, one has 

" ~ ^ 7 

1 — 

where (3o is the first coefficient of the /^-function (so /? = —(3oa‘l + ■ ■ with the obvious 
generalization to higher orders. Eq. (2.23) implies that 

[as,M]^/3oal (2.24) 


(2.23) 


Therefore, the explicit form of a general function of dg and M will depend on the choice 
of operator ordering. The ^-evolution equation eq. (2.5) becomes [24] 

doe. M) = x(a„ M)G(e. Af). (2.25) 

where the derivative with respect to M acts on everything to the right, and x may be 
expanded as in eq. (2.8) keeping the powers of dg on the left. It is easy to show that 
different choices of operator ordering in eq. (2.25) would corresponds to different choices 
for the argument of the running coupling [24,11]. 

When X is approximated by asXo? eq. (2.25) becomes, after taking a second Mellin 
transform, 

(1 - Poas:^)NG{N, M) + F{M) = asXo{M)G{N, M), (2.26) 

dm 

where the f un ction F{M) is a boundary condition. By expressing the anomalous dimen¬ 
sion in terms of the exact solution to eq. (2.26), it is possible to prove [7] that corrections 
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eq. (2.22) to naive duality eq. (2.9) which generalize the NLO term eq. (2.20) may be 
determined perturbatively, thereby ensuring that duality continues to hold at the running 
coupling level to all perturbative orders. Because the dual GLAP equation enjoys the 
standard factorization properties, this also implies perturbative factorization of the run¬ 
ning coupling BFKL equation. Furthermore, the all-order sum of these running coupling 
corrections can be determined explicitly when xo{M) is replaced by its quadratic approxi¬ 
mation near its minimum at M = 1/2: xo(Af) ~ c-|- ^k{M — i)^, because in such case the 
inverse Mellin eq. (2.6) of the solution to eq. (2.26) can be determined exactly in terms of 
Airy function of an argument that depends on as{t) and N. Since the asymptotic small x 
behaviour is determined by the central interval in M of x, one can show that the quadratic 
expansion leads to the correct asymptotic behaviour implied by the kernel xo- The all¬ 
order resummation of running coupling effects can then be obtained [7] by showing that 
this exact solution based on the quadratic approximation to x determines the asymptotic 
small X behaviour of the solution computed from the full x- 

The running coupling resummation of the DL anomalous dimension is found to greatly 
soften the asymptotic small x behaviour in comparison to the DL result with unresummed 
running coupling effects, by replacing the naive DL small x intercept with an effective 
one [10,7,8]. As a consequence, it is no longer necessary to fit the small x behaviour to the 
data in order to obtain reasonable phenomenology [7,8]. This, however, makes a precise 
description of the minimum of x mandatory. Specifically, at NLO of the DL expansion 
where no minimum is present there is a loss of predictivity because it must be assumed that 
the minimum is present in the all order result [8]. Fortunately, the minimum is restored 
at each order by exploiting the underlying BFKL symmetry, as we explain in the next 
section. 


3. Symmetrisation 

We have seen that for a precise treatment of non leading terms a better control of the 
central region near M = ^ is needed and for this one has to stabilize the bad behaviour of 
X near M = 1. This can be done using the symmetry of the BFKL kernel. The underlying 
Feynman diagrams which determine the kernel K ^cts, in eq. (2.7) order by order 
in perturbation theory are symmetric upon the exchange of the incoming and outgoing 
gluon. But the interchange of the incoming and outgoing gluon virtualities in the 

argument of the BFKL kernel K ^cts, in eq. (2.7) corresponds to the transformation 

M —1 — M in the argument of its Mellin transform M). Hence, x{(^s, M) must be 
symmetric upon this interchange. 

The symmetry, however, is broken by two effects: the running of the coupling and the 
fact that in deep-inelastic scattering the WMellin transform eq. (2.1) is performed with 
respect to the variable ^ = ln(s/(5^) which depends on the outgoing scale only. We will 
start with the fixed cts case, while running coupling effects and the associate symmetry 
breaking will be discussed in the next section. The kernel is then only asymmetric due 
to the choice of definition of As is well known [16], this asymmetry corresponds to a 
reshuffling of the arguments M, N of x^ namely, a symmetrization xe of the kernel x of 
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eq. (2.25) can be obtained from a symmetric kernel Xa which corresponds to the symmetric 
choice ^ = ln{s/^/Q^^) throngh the implicit eqnation [16]: 

XT,{as, M + \xa{M)) = Xa{as: M). (3.1) 

Hence, we can implement the symmetrization, thereby resnmming the M = 1 poles, 
by performing the donble-leading resnmmation of M = 0 poles of x, determining the 
associated Xcr throngh eq. (3.1), then symmetrizing it (as Xa mnst be symmetric), and 
hnally going back to DIS variables by nsing eq. (3.1) again in reverse. 

The procednre can be nnderstood most easily by hrst considering a toy case, where 
we pretend that the anomalons dimension exactly coincides with its leading-order form 
as 7 o(A^), so that the BFKL kernel is jnst given by Xs{ois/M) from eq. (2.18). This corre¬ 
sponds to a leading InQ^ (or collinear) approximation to the BFKL kernel. Symmetrization 
is enforced by nsing eq. (3.1) to relate the DIS kernel to the symmetric one: we replace M 
by M -|- A^/2 and symmetrize for M (1 — M). Note that this is allowed at the leading 
InQ^ level, becanse the symmetrizing terms are by constrnction free of poles at M = 0 
and are thns snbleading in InQ^. 

We thns obtain the following pair of kernels, related by eq. (3.1): 

Xatoy (as, M, N) ^ Xs( “w ) + 

\ ivi 2 2 ^ 

XT. toy (cts, M, N) — Xs (xf) + Xs ’ (3-3) 

where Xotoy is symmetric, and XTtoyicts, M) is relevant for DIS, i.e. dnal to the DIS 
anomalons dimension. 

Eqnations (3.2), (3.3) shonld be viewed as implicit dehnitions of the corresponding 
kernels Xo-(«s, -Tf), Xs(«s, M) as solntions of eqnations of the form 

x{as,M) = x{as,M,xia,,M)) (3.4) 

in respectively symmetric or asymmetric variables. So x{as, M, N) can be viewed as“off- 
shell” continnations of the kernels x(cts,-Tf), which are obtained from them by pntting 
N “on-shell” (see eq. (2.9)). Alternatively we can nse nse the off-shell kernels as implicit 
dehnitions of the anomalons dimensions 7 (as, N) throngh the dnality condition (2.17): 

N ^ x{(Xs,'y{(^s,N),N). (3.5) 

Snmmarising, we obtain x(cts, M) or 7 (as, N) by pntting either N or M “on-shell” (i.e. set 
N = x(cts, M) or M = 7 ( 0 : 5 , N) respectively) in the “off-shell” relation N = x(«s, Af, N). 

Note that a,s N ^ 00 , the symmetrising contribntion in eq. (3.3) Xs( i-m+at ) ~ 
Xs(^) ~ 0{^). Thns this term does not spoil the identihcation of the resnmmed kernel x 
with the GLAP resnlt Xs at large iV, which in tnrn gnarantees that the resnlting anomalons 
dimension matches smoothly to the GLAP anomalons dimension. 

However the symmetrizing contribntion does violate the moment um conservation 
condition Xstoy(cts,0) = 1, albeit by snbleading terms. Indeed, at the momentnm 



conservation point M = 0, Xs{c(s/M) by constrnction satisfies Xs{c(s/M)\m=o = 1- 
However, Xs{ois/{^ - M))|m=o = Xs(as) = 0{as), and Xs{ois/{^ - M + iV))|M=o = 
Xs(as)(l + 0{N)) = 0{as), where the last step follows from dnality eq. (2.9). Hence, the 
momentnm violation is 0 {as) = or next-to-leading order in an expansion 

of X in powers of cts at fixed asjM. A simple way of enforcing momentnm conservation 
withont spoiling the symmetrization properties is then to add to eqs. (3.2),(3.3) an extra 
term of the form 

Xmomias^N) ^ CmfmiN), (3.6) 

where /m(l) = 1 and /^(O) = /m(oo) = 0, for example 

■^""(^) ^ (ivTTF’ 

and with the constant Cm is chosen to fix np the momentnm constraint. The condition 
/m(0) = 0 ensnres that Xmom(«s,fV) is snbleading: becanse Cm = 0{as), Xmom(«s,fV) 
eq. (3.6) is of order 0{asN) [l + 0{as) + 0{N)]. Bnt by the dnality eq. (2.9) N = 
0{as/M), hence Xmom(cts,-^) is 0{al/M)^ i.e. next-to-leading InQ^. The condition 
limjv^oo fm{N) = 0 again ensnres that the extra term does not spoil the indentification of 
the resnmmed x kernel with the GLAP resnlt Xs at large N. 

The fnll symmetrized donble-leading resnlt can be obtained throngh a similar pro- 
cednre, bnt starting with Xdl LoiM), given in eq. (2.15). Clearly, the contribntion from 
the leading-order BFKL kernel xo is already symmetric, bnt the symmetry shonld only 
hold in symmetric variables, while in DIS variables it mnst be asymmetric. Fnrthermore 
the collinear and anti-collinear singnlarities mnst match exactly with those derived from 
GLAP. To achieve this it is necessary to modify the BFKL kernel by snbleading terms 
in the same way that we symmetrised the pnre GLAP model above. Recall that the LO 
BFKL kernel 

Xo(M) = ^ [2V^(1) - i;{M) - V^(l - M)], (3.8) 

where if;{M) = ■^lnF(a:) is the polygamma fnnction. We can separate it into a snm of 
two pieces: a collinear piece Xq^ with poles at M = 0, —1, —2 ,... coming from the collinear 
region > k‘^ and an anticollinear piece Xq with poles at M = 1, 2, 3,... coming from 
the anti-collinear region < k^: Xo — Xo + Xo with 

X([(M) = — (V^(l) - V^(M)) 

^ (3.9) 

Xo(M)^Xo+(l-^)- 

It follows that in symmetric variables the off-shell extension of xo ninst be 

W(M, N) = ^ [V^(l) + ^{1 + N)- ^{M + N/2) + N/2 )], (3.10) 

so that in DIS variables 

W(M - f, iV) = ^ [^(1) + ^(1 + AT) - ^(M) + iV)]. (3.11) 
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The leading collinear and anticollinear poles 


Xo"(M,iV) 


ric 1 I 1 

TT \_M+N/2 l-M+N/2\ ’ 


(3,12) 


then coincide with those from GLAP [eq. (3.2) expanded to 0 (q;s)], while at the same time 
the only difference between (3.10) and (3.8) is in terms of 0{N) = 0{as) [nsing (2.9) with 
x(M, CKs) = Q;sXo(Af)]. Note that besides the shift in the argnment of '0(1 — M) we also 
added a term to Xo proportional to '0(1 + N) — '0(1): this is to ensnre that at large N in 
DIS variables the anticollinear terms '0(1 + N) — '0(1 — M + N) ~ 0{1/N), and thns that 
the matching to GLAP in the collinear region (and in particular at large N) is not spoiled. 

Putting together the GLAP and BFKL components, the symmetrized off-shell leading 
order DL y is then 


XaLoicts.M, N) = X, 


M+ 


N 


+ Xs 


1-M+ 


N 


+ a,xoiM, N) + x mom 7 iV), (3.13) 


where if 

Xo(Af. N) = xJ(M. N) + XoiM, N), 

xJ(M.iV) = !;^(V>(l)-V'(AA+f)). (3.14) 

X^{M, N) = 2^(V>(1 + N)-i,{l-M+ f)). 

then 


Xo(M,iV) = xo(M,iV)-x^-(M,iV) 

(3.15) 

= ^ (0(1) + 0(1 + iV) - 0(1 + M + f) - 0(2 - M + f)) , 

where in the last step we have used the identity 0(a: -t- 1) = 0(a:) -|- i. In DIS variables we 
get 

= Xs (^) +Xs +asXo{M - f, iV) + Xmom(as, iV), (3.16) 

where 


Xo(M - f, iV) = ^ (0(1) + 0(1 + N)- 0(1 + M) - 0(2 - M + N)) . (3.17) 

This differs from the original unsymmetrized leading order DL result eq. (2.15) by sub¬ 
leading terms, as we now show. The symmetrizing contribution Xs[cts/(1 ~ Af)] is 0{as) 
and free of M = 0 poles, hence next-to-leading InQ^. Because it is an 0{as) contribution 
to X if is leading In A However, its 0{as) term is subtracted by the double counting 
term — ^0^ jir/gqrjY, and what remains is 0{a^g)^ namely next-to-leading In A Finally, the 
leading In - contribution to Xelo eq. (3.16) differs from xo because of the A^-dependent 
shift of the argument M: however, by duality eq. (2.9) N = 0{as) and thus the difference 
is 0{al)^ namely next-to-leading In-. The violation of momentum conservation in the 
absence of Xmom is still sub leading. This is proven by simply observing that Xmom(«s, N) 
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differs from the DL result by subleading terms. But the DL result respects momentum 
conservation up to subleading terms, and the remaining terms are subleading with respect 
to it everywhere. 

In a similar way, starting from xdlNL oiM, as) in eq. (2.16) , one can write down the 
symmetrized off-shell NLO DL y: 


XoNLoiois, M,N) — Xss i —) +Xss ( --^) +Xmom(cts, N). (3.18) 

\M+-^ / J 

In DIS variables this corresponds to: 

XJ:NLo{as,M,N) = (if)+Xss ( i_M+jv )+^^Xi {M - N)+Xmom{as, N). (3.19) 

The function xi which is related to the Fadin-Lipatov kernel xi (the next-to- 

leading contribution to the BFKL kernel eq. (2.8)) minus double counting, up to sub¬ 
leading corrections due to the symmetrization and to the running of the coupling, will be 
constructed explicitly in section 6 below. 

Let us now turn to the generic properties of the symmetrized on-shell kernels, obtained 
by solving eq. (3.4) in either symmetric or DIS variables. First, it is clear that, because 
moment um conservation in asymmetric (DIS) variables hxes x(cts,0) = 1, in symmetric 
variables the associate symmetric kernel will satisfy the constraint 


X<T = X<T (tts, I) = 1- (3.20) 

When transforming back to DIS variables, this will lead to 

XE(as, 0) = Xs(as, 2) = 1. (3.21) 

Since in the central region M ~ 2, x ~ 0{as) and thus rather less than one, at least for 
small enough values of as, on the real axis of the M-plane the symmetrized kernel will 
have a minimum order by order in the (symmetrised) double leading expansion. Since the 
change of variables from symmetric to DIS variables is just a shift M —M — there will 
be a minimum is DIS variables too. This is an important generic result, since it implies 
necessarily that at hxed coupling the splitting function must behave as a power of x, order 
by order in resummed perturbation theory. 

Now consider the position and curvature around the minimum. While Xcr has the 
minimum at M = 2, the minimum of xs is displaced, but the value at the minimum and 
the curvature are the same. To see this, we rewrite the implicit equation eq. (3.1) as 


Xa{(^s,M-\x^{M)) ^X^{(^s,M), (3.22) 

where Xa{ois,M) — Xcr(cts,l — M) is stationary at M = |. Eq. (3.22) implies that the 
derivatives of xe and Xo- are related by 


Xe(«s,M) 


Xe(«.,^) 


Xa'{as,M - |xe(M)) 

1 + iXa(«s, M - \x^{M)) ’ 
Xa"{as,M - |xe(M)) 

(1 + - ixE(M)))3’ 


(3.23) 
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Figure 1: Plot of different symmetric small x kernels x- From the top (in the central 
region): the LO BFKL kernel a^Xo, the NLO and LO resummed DL expansion kernels 
Xu, both on-shell and expressed in symmetric variables, the NLO BFKL kernel asXo + 
QsXi- All curves are determined with /3o = 0 (fixed coupling), as = 0.2 and u/ = 0. 
Note the unphysical branches of the LO and NLO BFKL curves, outside the interval 
0 < M < 1. 


where the prime denotes differentiation with respect to M. It follows that the minimnm 
of Xe(«s, Af) is at Mg, determined by 




(3.24) 


which implies 



(3.25) 


i.e., the intercept and cnrvatnre of xs and Xa at the respective minima are the same (note 
that this is not trne for higher derivatives of either function: after all, xe is asymmetric). 
Hence, whereas the location Mg of the minimum is shifted away from Mg = | by the 
asymmetric scale choice, the intercept and curvature at the minimum are unaffected. 

Although xdl has poles at M = ±1, ±2,..., in symmetric variables Xa is actually 
an entire function of M. To see this, note that the off-shell X(t{<^s, M, N) is of the form 
ip{M -t- y) + (/?(! — M + where (p{M) has poles on the negative real axis at M = —n, 
n = 1,2,..., so Xcr(cts, M, N) will have poles at M + ^ — —n, near which it must go to 
inhnity. Eq. (3.4) then shows that N ^ oo there, i.e., as we go on-shell the corresponding 
singularity in Xcr{oig,M) is shifted to the point at inhnity. Similarly, in DIS variables 
Xs(cts, Af, N) is of the form + (p{l — M -\- N), so by a similar argument the poles on 
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Figure 2: The resummed kernels of fig. 1 expressed in asymmetric variables as relevant 
for the DIS structure function evolution. Also shown for comparison the LO and NLO 
GLAP kernels and the LO and NLO double-leading (DL) kernels. 


the positive real axis are shifted to infinity when we go on-shell, and M) is free of 

singnlarities to the right of ReM = — 1. The singnlarities at M = — 1, —2,... correspond 
of conrse to higher twists, which are ontside the control of our approximations, but die off 
rapidly at high Q^. 

Note that whereas is an entire function, its perturbative expansion still has 

order by order singularities at M = 0, just like Xdl(«s, M) and Xe(«s, M). In fact whereas 
the perturbative expansions Xdl{c(s: M) and Xs(«s: M) have single collinear poles (of the 
form (cts/M)"^) the perturbative expansion of Xaicts^M) has double singularities of the 
form To see this substitute the expansion eq. (2.8) in eq. (3.1) to get 

Xi:o{M) =Xao{M), xei(M) = Xai(M) - ixo(M)xo(M),..., (3.26) 

whence the result. Similarly the anticollinear singularities in DIS variables, i.e in XE(cts, M) 
at M = 1 are also double singularities, of the form The spurious cubic 

singularity in xi is well known [16,25] and will be discussed in section 5 when extending 
the symmetrization to next-to-leading order. 

A comparison of the unresummed BFKL kernels ctsXo and ctsXo + ctsXi with the 
symmetrised DL expansions at LO and NLO is presented in fig. 1. The curves shown 
are, in symmetric variables, the on-shell versions obtained by using eq. (3.4), from 
eqs. (3.13), (3.18). First, we observe that the LO and NLO symmetrised DL curves are very 
close, indicating a good behaviour of the corresponding expansion, while the NLO BFKL 
kernel was in all respects completely different from its LO approximation. We also observe 
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that the DL expressions are very smooth and rather flat. This is dne to the fact that 
the symmetrization procednre shifts the poles at infinity while momentnm conservation 
imposes x{~\) =x(|) = 1 - fhe following section we will nse a qnadratic approxima¬ 
tion of the kernels near their minimnm. Clearly for the DL kernels the validity of this 
approximation is qnite accnrate in the whole physical region 0 < M < 1. The resnmma- 
tion considerably rednces the valne of x at the minimnm in comparison to the LO BFKL 
kernel, implying a softer small x asymptotic behavionr. Amnsingly, after resnmmation the 
NLO correction is actnally small and positive, in complete contrast to the large negative 
correction snggested by a naive interpretation of NLO BFKL. 

In fig. 2 we display a set of x kernels expressed in asymmetric variables as relevant 
for DIS strnctnre f un ction evolntion. The LO and NLO BFKL kernels, the LO and NLO 
DL cnrves and Anally those derived at LO and NLO from the symmetrization procednre 
are compared with the dnal of LO and NLO GLAP. The stabilization of the resnmmed 
expansion for M ^ ^ throngh the symmetrisation procednre is apparent: in fact it is 
gnaranteed by the anticollinear momentnm conservation constraint at M = 2. The position 
of the minimnm is shifted slightly to the right or M = i bnt its height and cnrvatnre are 
the same as in the resnmmed cnrves of fig. 1. However the resnmmed resnlt is now very 
close to the GLAP resnlt for all M below the minimnm of the resnmmed symmetrised 
cnrve. As we shall see in the next section, rnnning conpling effects extend this agreement 
considerably. 

Let ns Anally discnss the effect of the switch from symmetric to asymmetric variables 
on the expansion of the anomalons dimension. At the flxed-conpling level, this is simply 
determined nsing dnality eq (2.9) in eq. (3.22): if 7,7 and 7 e are respectively dnal to the 
symmetric and nonsymmetric Xcr and xe, then, to all pertnrbative orders 

XY.{as,N) ='^^{as,N)+ (3.27) 

This shows that the effect of the scale choice is snbleading, i.e. N~^). If is 

expanded according to eq. ( 2 . 12 ), then the term indnced by the scale choice (the last term 
in eq. (3.27)) is a contribntion to Xss of order {as/N)~^. Of conrse, if 7 is pertnrbative 
(i.e. if it vanishes in the cts — 0 limit), this term will cancel against an eqnal and opposite 

contribntion to x^ss- Fnrthermore, an important implication of eq. (3.27) is that the 
singnlarities of 7 ( 0 : 3 , A^) in the complex A^-plane are the same in both symmetric and 
asymmetric variables, and thns in particnlar the small-singnlarity which determines the 
small X behavionr of the splitting f un ction will also be the same. 


4. Running coupling 

So far, we have discussed the improvement of the x kernel and its dual anomalous 
dimension at the fixed coupling level. As discussed in the introduction, the main motivation 
for the symmetrization of the previous section is that the symmetrized double-leading result 
has a stable minimum order by order in the resnmmed expansion, in contrast to the BFKl 
or unsymmetrized DL expansions which have an instability. The existence of a minimum 
enables a running coupling resnmmation, which is required [7,8] if we wish to get a stable 
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resummed anomalous dimension at small N, and thus a uniformly stable splitting function 
at small x. To this purpose, in line with the strategy devised in refs. [7,8], the running 
coupling evolution equation (2.25) is solved with a quadratic approximation to the kernel 
near the minimum, but with the running of the coupling treated exactly at the leading 
\nQ^ level. The solution is then used to improve the anomalous dimension computed 
within a perturbative treatment of the running of the coupling. 

In refs. [7,8] the quadratic approximation near the minimum was applied to the LO 
BFKL kernel xo, eq. (3.8). This kernel is symmetric and has a minimum at M = i. 
However, we have seen that in order to preserve the symmetry and the existence of the 
minimum beyond the LO, we have to construct a different expansion for x which at LO is 
given by Xa Lo{ois^ M) obtained from its off-shell version eq. (3.16) through eq. (3.4). 

We are thus led to solve the equation 


7VG(7V, M) = xU^s, M)GiN, M) + F(M), (4.1) 


where Xai^siM) has been replaced by its quadratic approximation x^(ds,Tf) near the 
minimum 

xli^s: M) ^ c{as) + \K{as) {M - (4-2) 

and F{M) is a boundary condition to perturbative evolution. Note that, in practice, 
because of eq. (3.25), the intercept c{as) and curvature K{as) can be computed using 
the symmetric variable resummation Xcr- The quality of the quadratic approximation, 
especially in the relevant region 0 < M < 4^ can be appreciated from hg. 3, where 
the resummed curves obtained from the symmetrisation procedure, shown in hg. 1, are 
compared to their quadratic approximations. 

There are two complicating features in the present case with respect to our previous 
work [7]. First, in earlier papers the curvature c and intercept k of the kernel were linear in 
as', eq. (4.1) could then be solved to yield a solution expressed in terms of Airy functions. 
Here we need to expand c(ds) and ^(ds) in powers of dg — as and then, to the required 
accuracy, it turns out that the solution can be expressed in terms of Bateman functions. 
The second complication is that, since dg is an operator which does not commute with 
M, any expression like x(®s, M) is meaningful only after specifying the operator ordering. 
This amounts to specifying the argument of the running coupling. The hrst issue will 
be dealt with immediately, while the ordering issue will be postponed to the next two 
sections, where we will explicitly construct the resummed kernels, anomalous dimensions 
and splitting functions in the LO and NLO approximations. 

When, as in the present case, c(dg) and «:(ds) are given as power series in dg, the 
differential equation (4.1) cannot be solved in general. However, we can determine the 
solution by expanding dg about ctg to hrst order in dg/9o- It is easy to see that this 
procedure guarantees that the whole sequence of leading log contributions to Og is correctly 
included, up to subleading corrections. We have 


c(dg) 

K(dg) 


c(ag) -I- (dg 

K{as) + (dg 


- ctg) c'(ag) O 

- as) K\as) + O 



(4.3) 

(4.4) 
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Figure 3: Plot of the NLO and LO resummed DL expansion kernels Xcr, both on-shell 
and expressed in symmetric variables, compared with their quadratic approximations 
(dashed curves) near the minimum at M = i 


and, from eq. (4.1), 


N-c{a,) - ^K(a,) (M - i)^l G(N,M) = a, (c'(a,) + (M - i)h G{N,M)+F{M), 

(4.5) 

where we have dehned 


c{as) = c{as) - asc'{as) 

R{as) = K{as) - asK{as). 

Equation (4.5) can be simplihed dehning G{N, M) through the implicit equation 

N 


G{M, N) = 


N - c{as) - |k(q!s) (M - I) 


rG'(M, N), 


so G{N,M) satishes 


with the kernel 


NG{N, M) = as(f){M, N)G{N, M) + F(M), 
N{c'{as) + ^K'(as) (M - 


(4.6) 

(4.7) 


(4.8) 


¥>(VA/) = 


N — (e(as) + (M — 5 )' 


(4.9) 


(4.10) 
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Equation (4.10) has the form of a leading-order M, N space running coupling BFKL 
equation, generalized to the case in which the kernel (p{M, N) depends on both M and 
and not just on M as the BFKL kernel xo does. As discussed in ref. [8], the anomalous di¬ 
mension is entirely determined by the dependence of the inverse M-Mellin transform of 
the inhomogeneous solution to this equation, viewed as a differential equation in M. Fur¬ 
thermore, this inhomogeneous solution to all perturbative orders satishes a factorization 
property, whereby it can be written as the solution to the associate homogenous equa¬ 
tion, times an M-independent, but A^-dependent boundary condition. The anomalous 
dimension is therefore entirely determined by the solution to the associated homogeneous 
equation. The fact that the kernel now depends also on N does not affect the form of the 
solution, since the inhomogeneous equation is a differential equation in M. Hence, using 
the general solution of ref. [7] (given as eq. (3.21) of that reference) we get 


G{N,t) - 


dM 


^—ioo StTZ 


exp 


M 


r-M 


(3oas(t) (3o 


- I dM'<P{M\N) 


^ gl/2/3oaAi) 




(4,11) 


27ri \l — A{asTN)m 




where we have dehned 


while 


A(a.,iV) 


S(a„iV) 




(N-c(as))’ 


1 / c'(a^) 

2/^0 \JV-c(a^) 


1 _ 1 

as(t) as(t) 




(4.12) 

K'(as)\ 

In - c{ag) 

(4.13) 

K{ag) ) 

V ’ 

K'{ag) 


(4.14) 

K{ag) ' 



The Mellin inversion integral can now be computed exactly, with the result 


G(N, t) = (3oas{t)B{as, N)K2 b{c,,n) {1/Po(^s{t)A{as, N)) (4.15) 

where K^j{x) is the Bateman function, dehned as the solution of the differential equation 

-A''(a))+(l-^)A, = 0, (4.16) 

with boundary condition Ar,^(0) = 1. Note that the argument of A and B eqs. (4.12)- 
(4.13) is the hxed coupling cts, not agit). This is due to the fact that in the approximation 
eq. (4.3),(4.4), all functions of ag are linearized by expanding about the hxed coupling ag. 

The anomalous dimension is determined by taking the logarithmic derivative of 
G{N,t), and then setting as{t) — ag. The extra factor which relates G to G eq. (4.8) 
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only affects the boundary condition which, as proven in ref. [7], does not contribute to the 
anomalous dimension. Therefore, we get 


7b (as, iV) = ^ In G{N, t) 


di)=c 


— k ~ PqO^s + 


K 


2B(as,Ar)(^//^oas^(as, N)) 
A{asjN) K2B(a^,N)i^/Po^sA{asj N)) 


(4.17) 


The Bateman anomalous dimension can be matched to the naive (i.e. hxed coupling) 
dual of the x with the same procedure which was used in ref. [7] for the Airy anomalous 
dimension. Namely, we determine 7 b( as, N) in the /3o —^ 0 limit by saddle-point evaluation 
of the Mellin integral eq. (4.11). The leading order saddle condition is 


7V = a,(t)(p(iV,Ms), 
where <f>{N,M) is given by eq. (4.10), which reduces to 

N = Cs{as{t), as) + |Ks(as(t), as)(Ms - ^)^, 

where 

Cs(as(t), as) = as(t)c'(as) + c(as) 
Ks{as{t), as) = as{t)K'{as) + it(as) 


(4.18) 


(4.19) 


(4.20) 

(4.21) 


are the linearized forms of c and k. Solving eq. (4.19) for Ms{N) gives 


Ms(iV) = i- 


( iV - c(as(t),as)) 
i«;(as(t), as) 


which, if one lets as(t) = as, reduces to the standard naive dual 


7f (as,iV) = ^- 


{N - c(as)) 
hi^iois) 


(4.22) 


(4.23) 


of Xai^sjN) eq. (4.2) as it should. 

Further terms in the expansion about the saddle generate an expansion of the anoma¬ 
lous dimension in powers of particular, the fluctuations about the saddle 

lead to a contribution to the anomalous dimension which can be computed as in ref. [7]: 


1 

(iV,Ms(as(t),as,iV)) 


exp 


dt' Axss{oisit'), as, N) 


'to 


1 


(iV,Ms(as(t),as,iV))’ 

(4.24) 
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Figure 4: The Bateman anomalous dimension 'yBiois, N), computed with as = 0.2 and 
the values of the parameters which correspond to the LO resummed curve of fig. 1, 
namely c = 0.27, c' = 0.71, c = 0.13, k = 1.40, k' = 0.99, R, = 1.20. 


where the prime denotes differentiation with respect to M (which is then set eqnal to Mg, 
eq. (4.23)), and 


{as{t),as,N) 


-\(3oas{t) 


M) 


M=Ms 


(4.25) 


Use of the saddle condition eq. (4.22) with the kernel eq. (4.10) in eq. (4.24) gives, after 
setting as{t) = cts, 


7,1 (a., N] = (a., N) + (4,26) 

where 

7^,0 N) = -Poag + (4.27) 

Note that 7 ^ depends on both ag/N and agi in fact, it differs by terms of order 0{as) 
from a contribntion of order 'jgg in the expansion eq. ( 2 . 12 ) of 75 in powers of ag/N at fixed 
N. The contribntion 7 ^ 0 ( 0 )^,^) = limjv^oo 7 ^ (cts, A^) is snbleading in the expansion 
eq. (2.12) bnt leading-order in the expansion eq. (2.3) of the anomalons dimension in 
powers of ag at fixed N. It follows that, at LO in the symmetrised DL expansion, as 
N ^ 00 

7B(a„ N) ~ 7 f (a., N) + 7f,,o(a.. V) + 0( 1/Af). (4.28) 

The Bateman anomalons dimension eq. (4.17) is displayed in fig. 4, with the valnes 
of the parameters which correspond to the qnadratic approximation to the LO resnmmed 
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kernel displayed in figs. 1,3. It is apparent that its behavionr is similar to that of the 
Airy anomalons dimension discnssed in ref. [7]: the anomalons dimension is regnlar and 
monotonically decreasing along the real axis, with a simple pole at = Nb{c(s), located at 
the rightmost zero of the Bateman fnnction Kj^{x). The small x behavionr is determined by 
the location of this pole, which in tnrn depends on the valnes of the intercept and cnrvatnre 
of the qnadratic kernel c{as) and K{as) eq. (4.2) and their hrst derivatives c'(as), K'{as)- 

The qnalitative behavionr of the Bateman rnnning conpling resnmmation can be in¬ 
ferred from £g. 5, where we compare the qnadratic approximation to the LO resnmmed 
kernel (same as in hg. 3) with the naive dnal eq. (2.9) of the physical branch of the Bate¬ 
man anomalons dimension (i.e. with N > Nb{c(s))- For comparison, the corresponding 
Airy anomalons dimension is also shown. The latter is only dehned when c and k are 
linear fnnctions of as, and it is thns determined by taking in eq. (4.2)c(ds) ~ dsc(as)/as 
and K{as) ~ ds«:(as)/as. This is a worse approximation than that nsed to determined the 
Bateman anomalons dimension, since the leading log terms are not correctly reprodnced. 
The resnmmation of rnnning conpling effects effectively takes taking the minimnm of x 
from one half ont to inhnity [7], thereby greatly extending the range of the physical region. 
For the anomalons dimension this means that the nsnal hxed conpling cnt at = c(as) is 
replaced by a simple pole at A^B(as) to the left, i.e. NB{ois) < c(as). This in tnrn means 
that the small x behavionr of the splitting f un ction will be softer (more slowly rising) than 
that expected from naive dnality. 

In hg. 5 (right) we display the Bateman effective kernel determined after snbtracting 
from the Bateman anomalons dimension its A^ —cx) limit, namely, the dnal of 75 — 7 ® q, 
with 7 ^ Q given by eq. (4.27), together with the corresponding Airy curve. This comparison 
shows that, after this constant subtraction, for M < 0 7 ^ is very well approximated by 
the leading order asymptotic /I —0 approximation eq. (4.23), i.e., by the fixed-coupling 
quadratic curve, as is 7 ^ 1 . In other words, the running coupling resnmmation for N > c(as) 
(to the right of the hxed-coupling cut) is a very small correction. 

The location of the rightmost singularity of the anomalous dimension, which deter¬ 
mines the small x behaviour, is shown in hg. 6 as a function of ag. This singularity is a 
cut at hxed coupling, and a simple pole after running coupling resnmmation. The nonlin¬ 
ear dependence of c on a^ in the resummed hxed-coupling case (as opposed to the linear 
dependence in the unresummed BFKL LO case) is apparent, and it explains why a linear 
extrapolation from the origin would not be a good approximation for realistic values of as- 

In all cases, running coupling resnmmation softens the leading singularity. The ehect 
is most dramatic for the LO BFKL kernel dgXo eq. (3.8), whose running coupling resum¬ 
mation is ehected by the Airy method of ref. [7]. Even though the percentage softening 
is more moderate for the Bateman resnmmation discussed here (see also hg. 5), the sym- 
metrization discussed in the previous section also reduces the hxed coupling intercept in 
comparison to the BFKL one. Hence, after running coupling resnmmation the LO linear 
(Airy) result and the NLO resummed (Bateman) result for the location of the pole end up 
being quite close. The stability of the resummed result when going from LO to NLO is 
further enhanced by running coupling resnmmation, so all three running coupling curves 
end up being close to each other. Interestingly, the contributions to the NLO kernel pro¬ 
portional to (3o (to be discussed in detail in section 6 ) further reduce the diherence between 
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Figure 5: (left) The quadratic kernel (upper curve) obtained from the expansion near the 
minimum of the LO resummed curve in fig. 3, compared with the corresponding effective 
Bateman kernel. The effective Airy kernel, also shown for comparison, is determined 
with the same valne of c and k, i.e. approximating c{a.s) ~ 6isc{as) / ots and ^(ds) ~ 
asK,{as)/as. 

(right) Same as in the left fignre, but with the effective kernels determined from the 
Bateman and Airy anomalons dimension after subtraction of the asymptotic large N 
contribution 


the resummed LO and NLO minima, and in fact they pull the NLO resul below the LO one 
for as ^ 0.15. For comparison, we also show the location of the NLO unresummed BFKL 
stationary point, though this is in fact unphysical for realistic ag ~ 0.1 where the NLO 
BFKL kernel in fact has a minimum instead of a maximum, so the associated anomalous 
dimension becomes imaginary and the cross section oscillates [26,27]. 

The Bateman anomalous dimension will be used in the next sections to construct the 
improved anomalous dimensions in the symmetrised double leading expansion, hrst at LO 
and then at NLO. 


5. Leading Order Resummation 

We wish now to determine anomalous dimensions and splitting functions by combin¬ 
ing the ingredients of the previous two sections: symmetrization of the double-leading 
expansion and its running coupling resummation. In previous work [7,6] we constructed 
resummed anomalous dimensions by exploiting the fact that the DL expansions of y and 
7 are dual of each other, up to subleading terms. We thus constructed the anomalous 
dimension 7 starting from its own DL expansion, namely 

iDhiots, N) — 'jDLLoiois^ N) -b as'jDL NLo{ois, N) + . . . 
iDLLoi^s, N) = [as'yoiN) + 7s(^) - wi (5.1) 

/ ^2 61 \ 

lDLNLo{as,N) = [as7i(iV)+7s^(^)-eo - q ;^ VTW ^ iVW 
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Figure 6: The leading (rightmost) singularity of various anomalous dimension as a 
function of as, all determined with n/ = 0. BFKL LO , res fix LO and res fix NLO 
(dashed) denote respectively the location of the cnt for asXo eq. (3.8) and for the f3o = 0 
resnmmed kernels x<t,lo and x<t,nlo obtained pntting on shell the kernels eq. (3.16) 
and (3.19). LO Airy, LO Bateman and NLO Bateman (solid curves) denote the location 
of the pole Nb for the running coupling resummation of the corresponding three curves. 
BFKL NLO (dashed) denotes the location of the stationary point of the NLO kernel 
eq. (2.8). 


Here, instead, in order to take advantage of the improved kernel which is 

the on-shell version of eq. (3.16), we mnst dehne the anomalons dimension Xt, lo{c(s, N) 
nsing the dnality eqnation (3.5): 

XT,Lo{ois,'y^Lo{ois,N)) = N. (5.2) 

This coincides with the LO donble-leading anomalons dimension Tdllo up to snbleading 
terms. 

To go beyond the hxed conpling limit, in the kernel xt,{c(s, M) we mnst make the 
replacement ag —> dg. As proven in Ref. [7], factorized dnality remains valid order by 
order in pertnrbation theory at the running coupling level, i.e. for any given y there exists 
a 7 such that the solution to eq. (2.26) with the given y and the solution to the usual 
evolution equation ( 2 . 2 ) coincide if, order by order, the boundary conditions are suitably 
matched and appropriate /3o-dependent corrections are added to the hxed coupling kernels. 
In particular, in the LO expression of xt.{c(s, M) eq. (3.16), the terms originating from 
Xo are not modihed and, given the basic commutator relation, eq. (2.24), they do not 
depend on operator ordering up to NLO corrections. For the Xs terms it is easy to see 
that an operator ordering exists such that the Xs which is dual to 70 coincides with the 
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fixed-coupling expression eq. (2.18). Indeed, the leading-order form of the duality relation 
eq. (2.17) in the running coupling case is 

7o(Xs) = (5.3) 

which immediately implies that if Xs satishes eq. (2.18), i.e. it is the inverse function of 
7o, then 

Xs [70 ((dj'M)-!)] = Xs = N (5.4) 

by direct substitution of eq. (5.3). Hence, the running-coupling dual of 70 coincides with the 
hxed coupling dual Xs, provided Xs is now evaluated as a function of the operator M~^as- 
A more detailed discussion of operator ordering issues in running coupling duality will be 
given elswehere [28]. 

Let us now turn to the symmetrization of Xs- First, we should observe that the 
symmetry of the underlying Feynman diagrams which determines the symmetry of the 

BFKL kernel K ^cts, j eq. (2.7) holds for the unintegrated gluon distribution Q^), 
obtained by differentiation from the standard gluon distribution which enters the GLAP 
equation: 

= (5,5) 

i.e., upon Mellin transformation, 

^(e,M) = MG(e,M). (5.6) 

When the coupling runs, the M-space evolution kernel acquires a further contribution when 
switching from the integrated to unintegrated gluon distribution due to the commutator 
of ds and the M prefactor in eq. (5.6). Indeed, 

X. {M-^as) M) = Xs {M-^as) M) 

So, going from the integrated to the unintegrated distribution switches the ordering of dg 
and M~^ in the argument of Xs- Note that this transformation is manifestly not symmetric 

upon M —> 1 — M, hence if the kernel at the unintegrated level is symmetric, the kernel 

at the integrated level won’t be. 

At the running coupling level, the kernel is symmetric only if the argument of the 
running coupling is also treated symmetrically. In fact, upon Mellin transformation, a ^ 
evolution equation of the form 

becomes 

^G(e, M) = [aJx?’(Af) + xt^{M)ai\ G(i, M), (5.9) 

p=l 
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where Xr\m) are the Mellin transforms eq. (2.6) of the kernels K^\q‘^/ k"^) 

and /k"^)^ respectively. Becanse the symmetrization is based on the interchange 

Q'^ /c^, it follows that the symmetrized version of a’lf^M) is 2 (^agf{M) + /(I — M)a^Y 
Fnrthermore, becanse 

= —PoM~^alM~^ = —(3oasM~‘^asj (5.10) 

while 

[d„ (1 - M)-1] = Po{l - - M)-1 = -/lo«.(l - M)-2d„ (5.11) 

it is easy to show that if 

/(d,M-i) = fiM-^as) = J2fk^sM-^ 

k 

then 

/((i - M)-'d,) = /(d,(i - M)-‘) = Ml - mM&I 

k 

It follows that at the rnnning-conpling level the symmetrized LO kernel is 

XaLo(&s,M,N) = (ds (M + +Xs ((l - M + .f-) ' Q',)+Xo(d,. Af. Af)+Xmom, 

(5.12) 

where 

Xo{as,M,N) = as(xQ{M,N) - —^) + (xq{M,N) - _ J 

= ^ [d. (</>(!) - </>(! + M + f)) + (V>(1 + JV) - V'(2 - M + f)) dj , 

(5.13) 

which are the rnnning conpling analognes of eqs. (3.13), (3.14). Different orderings of the 
Xs term can be determined nsing the commntators (5.10) and (5.11). Different orderings 
of the remaining terms are eqnivalent, np to snbleading corrections which modify the form 
of the snbleading xi, and will be discnssed in the next section. Here, we have chosen the 
operator ordering in the xo terms snch that it matches the ordering of the corresponding 
Xs terms, namely, so that the ordering is the same in the donble connting contribntions. 
Note that this corresponds to taking, after inversion of the Mellin transform, cXsiQ'^) in 
the collinear donble-connting term, and as{k‘^) in its anti-collinear connterpart. 

Starting from eq. (5.12) one can expand Xa Loi^si N) in powers of (M — 2) at 
hxed N and dg, aronnd its minimnm at M = 2 with a resnlt of the form: 

X.(ds, M, N) = x^(d„ M, N) + xl(d., M, N) + 0(/?od.) + 0((M - i)^), (5.14) 

where 

xli&s.M.N) = as (co(iV) + \ko{N){M - , (5.15) 
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while 


(5.16) 


X^(d„ M, N) = c(a„ N) + lR{a,, iV) (M - i)^ 
and the expansion coefficients are given by 


c{a,,N) = 2c, (iflv) ’ 

^(ds, N) = 2 Kg ^ 1-1-tv) ’ 


(5.17) 


in terms of the coefficients of the Taylor expansions of Xs'- 

( m) ^ Cs(as) + 5s(as) (M - i) + i«;s(as) (M - + ... (5.18) 

The on-shell qnadratic expansion eq. (4.2) can now be obtained by pntting N on shell, i.e. 
by solving N = Xai^s, M, N) with the right-hand side given by eq. (5.14) for N consistently 
to order (M — 4)2^ then identifying the resnlt N = Xo-(ds,M) . Alternatively one of 
conrse conld first obtain the fnll on-shell Xcr(ds, M) by solving eq. (3.4) and then expanding 
the resnlt abont M = A the resnlt will be the same. 

The 0(/3ods) corrections in eq. (5.14) are dne to operator ordering, and will be dis- 
cnssed in sect. 6. They are of relative order a, with respect to the cnrvatnre and intercept, 
which are determined to all orders in a, in order to obtain a good description of the x 
kernel in the vicinity of the minimnm. Fnrther O{(3oas) corrections, also to be discnssed in 
the next section, are generated when switching from the nnintegrated distribntion (which 
evolves with a symmetric kernel) to the standard integrated one which enters the GLAP 
eqnation. 

The final expression for the improved anomalons dimension at LO is now obtained 
by implementing the rnnning conpling resnmmation throngh the qnadratic approximation. 
Namely, we write the evolntion eqnation eq. (4.1) as 


7VG(7V, M) = Xa(«s, M)G(iV, M) + F(M), (5.19) 


and then exploit the symmetry of Xo- to expand it abont its minimnm, np to a qnartic 
term which leads to asymptotically snbleading [7,8] corrections. The anomalons dimension 
determined from solntion of eq. (5.19) in the qnadratic approximation is then simply given 
by 7 b eq. (4.17), bnt with the replacement | ^ ^ . The anomalons dimension in 

asymmetric (DIS) variables is obtained from it nsing eq. (3.27), which simply restores M, 
to the valne eq. (4.17). 

Hence, the resnmmed anomalons dimension after rnnning conpling resnmmation is 
given by 

7e LO if )) = 7E LO [N, a, [t )) + 7 ^ (a^ (t), iV) 

^ „ (5-20) 

- 7f N) - xf{as{t), N) + 7mom(as(t), iV), 

where 7 elo (-^7 cts(f)) is determined from eq. (5.2), 7 ^(as(t), A^) is the Bateman anoma¬ 
lons dimension eqs. (4.17), and xf ic^sf), N), Xo i<^sit), N), obtained letting a, otsf) 
in eqs. (4.23),(4.27), remove the donble connting between the the first two terms. Notice 
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that 7^ removes the square-root branch cut of the hxed-coupling anomalous dimension, 
thanks to the fact that the curvature and intercept of the quadratic kernel are the same 
eq. (3.25) in symmetric variables (used to compute 75) and DIS variables (used to compute 
7x;). The momentum subtraction term is needed to cancel small (and formally subleading) 
momentum violations incurred during the running coupling resummation: it can be chosen 
as 

7mom(as(f), N) = dmfm{N), (5.21) 

where fm must satisfy /m(l) = 1 and /m(oo) = 0 [8] and can thus be implemented through 
the same function fm introduced in eq. (3.7), and dm is a constant of order ct^. 


6. Next-to-Leading Order Resummation 

The extension of our results to next-to-leading order requires two steps: hrst, the 
symmetrization of the next-to-leading DL kernel eq. (2.16), and second, the next-to- 
leading treatment of the running coupling resummation. The general structure of the 
next-to-leading order symmetrized kernel was already given in eqs. (3.18),(3.19) in the 
hxed coupling case. At the running coupling level, operator ordering in the leading order 
contribution affects the form of the NLO term. We assume therefore that (in symmetric 
variables) 

X^(ds,M, iV) = XaLo{dis,M,N) + a^Xa NLoi^s, M, N), (6.1) 

where X(tL o{d(s, M, N) is given by eq. (5.12). The corresponding kernel in DIS variables 
is 

XE(d., M, N) = M-^,N). (6.2) 

Note that at the running coupling level the kernel for evolution of the integrated and 
unintegrated parton distributions are not equal, as shown in eq. (5.7) for Xs, while in 
fact, as discussed in the previous section, only the unintegrated distribution is symmetric. 
Hence, we will construct the symmetrized x kernel at the unintegrated level, and then 
switch back to the integrated level at the end when determining the resummed GLAP 
anomalous dimension from it. 

At the leading order level, we have already shown in section 3 that the hxed coupling 
XELo(as, Af, N) =XaLo{c(s,M-^,N) is equal to xdl Lo(Af, N) up to subleading terms, 
and generalized this in section 4 to the running coupling xt, Loicng, M, N). We must now 
check this equality up to NLO at the running coupling level, and compute the correction 
terms which are necessary to restore it when violated. This will lead to the determination 
of Xt. NLo{dis, M, N) eq. (6.2), and in particular of the function xi(Af) which was left 
implicit in eqs. (3.18),(3.19). 

We start therefore with xe Lo(ds, Af, A^) given by eq. (5.12). This differs from the 
double leading xdl lo{M, N) in two respects: because of the addition of the anticollinear 
term Xs((l ~ Af -|- A^)“^ds), and because xo is put off shell. In order to compare to the 
DL result it is convenient to view the simple and double poles {1 — M + N)~^as and 
((1 — M + N)~^as)‘^ of the anticollinear term as contributions to the off-shell xo and xi, 
respectively. 
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Consider first the anticollinear term Xs((l — M + N)~^as) of eq. (5.12). This term 
starts at 0{as), but the hrst two terms are already included in xo and Xh hence the 
remaining contribution starts at 0(af), namely NNLO in an expansion in powers of as 
at hxed M, and N^LO in an expansion in powers of as at hxed as/M. Let us turn now 
to the rather more complicated case of the off-shell extension of xo{M), which is given by 
the 0(ds) contribution to xe Lo(ds, M, N) (in an expansion in powers of as at hxed M). 
It is thus 

asX^iM - f, iV) + XoiM - f, N)as (6.3) 

- d, [xUM - f, iV) + XoiM - f, N)] - al/3o^i;\l - M) + 0{al), (6.4) 

where x^{M,N) are given by eq. (3.14), and in the last step we used the commutator 
eq. (2.24), which implies 

[a.. f(M)] = alliof'(M) - + ... (6.5) 

The ordering of eq. (6.4) for the 0{as) contribution in the expansion eq. (2.8) of the kernel 
has been adopted in refs. [5-6] for the construction of the double-leading expansion, though 
the ordering eq. (6.3), as discussed in sect. 4, is more convenient in a symmetrized approach. 

Let us hrst consider the ehect of the oh-shell extension on terms which are next-to- 
leading in InQ^, i.e. in the expansion of x ia powers of as at hxed asjM. The ehect of 
the oh-shell extension can be determined by considering 

[x3‘(M) + Xo (M - Af)] - xo(M) = ^ |</.(1 - M) - </.(! - M + A^)] 

= V>(1 + Af)) - M |V.'(1) - V>'(1 + Af)]} + 0{M^), 

where x^iM) are given by eq. (3.9). It follows that 

asxtiM - f , iV) + XoiM - f , N)as - dsXo(M) = 0(M) + ^(d^), 

where now XoiM,N) is given by eq. (3.14), and the 0{a‘^) correction eq. (6.4) due to 
operator ordering will be included in xi as we discuss below. Therefore, thanks to the 
replacement '0(1) —> '0(1-1-A^) in x~iM, N) eq. (3.14) in comparison to x~iM) eq. (3.9), the 
combination eq. (6.4) which appears in the leading-order kernel Xelo, eq. (3.16), dihers 
from the BFKL kernel ctsXo, eq. (3.8), by terms which are of order 0{asM) = 0{a‘^{^)~^) 
i.e. next-to-next-to leading InQ^ in an expansion in powers of as at hxed as/M. 

Next we consider the ehect of the oh-shell extension on terms which are next-to-leading 
In x, i.e. the expansion of the kernel in powers of ag at hxed N. Because this expansion 
of the kernel (and in particular the NLO term xiiM)) is usually computed in symmetric 
variables, it is useful to compare the combination which appears in the leading-order kernel 
Xelo, eq. (3.16) with the LO BFKL kernel ctsXo in symmetric variables, too. It turns out 
that the oh-shell extension of xo induces next-to-leading In ^ contributions which must be 
subtracted. Indeed, 

[xJ(M, iV) +X0(M, iV)] -xo(M) = [20'(1) - 0'(M) - 0'(1 - M)] N + OiN^). (6.8) 


( 6 . 6 ) 


(6.7) 
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In order to construct the appropriate subtraction, start with the NLO BFKL kernel 
eq. (2.8), written in symmetric variables. The NLO version of eq. (6.3) is (in symmet¬ 
ric variables) 

XaNLo{as,M,N) = asXo{M,N) + Xo {M, N)as + ajxiiM, N) -i-O(d^), (6.9) 

where the ordering in the NLO term xi{M) is immaterial since it only affects NNLO terms. 

With the symmetric ordering of dg adopted in eq. (6.4), the NLO contribution 
Xi(M, N) is symmetric. The on-shell contribution xi{M) has been determined in refs. [16- 
19], in the form 

Xi&s, M) = asXo{M) + alxi^{M) + O(d^) (6.10) 

= d.x+(M) + XoiM)as + [xr(M) + - M)] + 0(d3),(6.11) 

2 

where Xi^iM) is given by ^(y) in ref. [16], with the identification 7 —M and Xi^ = 

In practice, the effect of the Po term in eq. (6.11) is thus to reverse the sign of the ppi — M) 
contribution to A( 7 ) (as defined in ref. [16]) thereby leading to a manifestly symmetric 
result for 

Xi(M) ^ x(HM) + - M). (6.12) 

We can now finally remove the terms generated by the off-shell extension eq. ( 6 . 8 ) up to 
NLO by defining (in asymmetric variables) 

XNLLx{as,M) = asxtiM) + Xo {M - N)as + d^Xi(M) +0(d^) (6.13) 

Xi(M) = xi(M) - ^Xo(M)^ [2P'{1) - P'{M) - ^'{1 - M )]. (6.14) 

With this definition, the off-shell extension eq. (6.13) coincides with the on-shell 0{a‘l) 
kernel of refs. [16-19], up to sub-subleading corrections. 

Note that although xi(Tf) has triple poles at M = 0 and M = 1, these are cancelled 

in xi(-Tf), so that the leading small M singularity is O (a/s) as required by the matching 
to Xs- This cancellation would be spoiled if one did not replace xi with xi on the r.h.s. 
of eq. (6.14). Of course this is all as it should be: we obtained xi through eq. ( 6 . 8 ) as 
compensation for the off-shell extension of xo eq. (3.14), which was in itself constructed 
precisely so that LO collinear evolution is not spoiled at large N. The cancellation of the 
triple poles is then the NLLx manifestation of this fact. 

Therefore, to match to the collinear contributions Xs and Xss, and their anticollinear 
counterparts, as described by eq. (3.13) and eq. (3.18), we define a function x{M) as 


Xi(M) = Xi(M) 


I M 1-M 


+ 92 



(1-M)2 


(6.15) 


with coefficients gi are determined by demanding regularity as M —0 and M —1, and 
thus correct matching to Xs and Xss'- 


91 

92 


rifric 



lln^ 3nf 

127r2 TFTT' 


(6.16) 
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Now, note that, of course, a shift by an amount proportional to N of the argument 
of xi is sub-subleading: Xi(M -|- kN) = xi(M) -|- 0{as)- Hence, we can simply take in 
eqs. (3.18),(3.19) 


Xi(M, N) = Xi(M - f, N) + Oias) - Xi(M), (6.17) 

with xi(-^) given by eq. (6.15). With this choice, obtained putting on shell 

XT^iM, N) coincides with its double-leading counterpart up to subleading terms: eq. (6.7) 
shows that the off-shell extension of xo only adds terms which are subleading in an expan¬ 
sion in powers of cts at hxed cts/M, and eq. (6.14) shows that it only adds terms which are 
subleading in an expansion in powers of cts at hxed M, while the anticollinear contribution 
is also subleading. 

The choice of xi eq. (6.17) has two shortcomings, however. First, eq. (6.17) implies 
that the pair of kernels in symmetric and asymmetric variables are only related by the 
shift eq. (3.22) up to subleading corrections as well, because the contribution xi does 
not satisfy eq. (3.22). This in particular means that also the equality of curvature and 
intercept eq. (3.25) of the pair of kernels is spoiled by subleading terms. As discussed in 
sect. 4, it then follows that the cancellation of the square-root cut in eq. (5.20) is spoiled 
by subleading corrections. Second, as discussed in sect. 3, xe in symmetric variables is an 
entire function, provided that eq. (3.22) holds. This property is lost if eq. (6.17) is used, 
and Xe then has subleading poles on the real axis, so that in particular the quadratic 
approximation is much poorer. 

It is therefore necessary to construct xf{M,N) such that while it is symmetric in 
symmetric variables, in asymmetric variables xti^ ~ reproduces the correct LO 

and NLO singularities, as in eq. (6.15), even at large N, while at small N xi{M,N) = 
Xi{M) + 0{N). The construction follows the same procedure as at LO, in particular by 
hrst separating 

Xi(M)=.x^(M) + xrW, (6.18) 

where Xi"(-^7 -^) N)) is regular when ReM > ^ ( ReM < ^) just as in eq. (3.14). 

However the details are somewhat technical, and have been relegated to an appendix. 

Because operator ordering issues are only relevant in the leading order term,^ this 
concludes the construction of xe eq. (6.2): the NLO kernel Xt,nlo is found using the kernel 
Xi{M, N), which reduces to xi(Af) as 0 in eqs. (3.18) or (3.19), while the collinear 
and anticollinear terms are given by Xss eq. (2.19), with argument M~^as in the collinear 
and asM~^ in the anticollinear terms, and off-shell extension as in eqs. (3.18),(3.19). 

Let us now turn to the running-coupling resummation. First, we must determine the 
quadratic approximation to the NLO kernel obtained by putting on shell xe eq. (6.1). 
Operator ordering issues are only relevant in the leading-order terms, as otherwise they 

^ Strictly speaking, at the running coupling level the expression eq. (2.19) of Xss is corrected 
by 0(/3o) terms which are of the same order as Xss (5.10). However, since our final aim is the 
determination of the anomalous dimension, we will neglect these terms, which do not affect the 
final double-leading anomalous dimension, but only the NLO symmetrization. 
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lead to sub-subleading corrections, whereas different choices of operator ordering in xe lo 
lead to terms of order Xt,nlo- 

The leading-order kernel Xa lo is determined with the canonical ordering eq. (5.12), 
whereas the quadratic kernel eq. (4.2) used in the running-coupling BFKL equation (4.5) 
is based on the ordering eq. (6.4), with all a a to the left. The quadratic approximation 
thus receives two contributions due to operator reordering. The first is due to reordering 
in xoiMjN) eq. (5.13), analogous to eq. (6.4): 

M, N) = a.MM, N) + N) + O(filal) 

xt(M,N) = - Af + f) - ; N -M + f). 


where Xo{M, N) is given by eq. (3.14). 

The second is due to ordering in the collinear and anticollinear contributions Xs in 
eq. (5.12), which can be determined by expanding about M = 4, just as in eq. (5.14), but 
now retaining commutator contributions up to order dg. Note that, because the symmetry 
is broken by operator ordering, we now also get a cubic contribution. We obtain 

Xs {aa{M + f )“^) + Xs{il- M + ^)~^)aa = xli^s, M, N) xf° (a^, M,N)+0 ((ds/3o)^) 

( 6 . 20 ) 

where xti^s, M, N) is given by eq. (5.16) and the commutator correction is 


x1°(&.,M,N) = c^’‘(a„N)+P’‘{&„N) (M - + (M - l)"+0((A/-i)='), 

( 6 . 21 ) 

where 


c^°(«s,iV) = -2/lo: 
j^°(ds,iV) = -/3o 


a 


^ c' 


ac 


1 + iv \i + N r 


a. 




k^°(ds,iV)--2/lo- 


(l + iV)2 \1 + Nj^ 
2 


a. 


( 6 . 22 ) 


ct: 


8 c' 


(l + iV)3 [ \1 + N 


etc 


46' 


Ots \ / 


etc 


14-iV 


where the functions c(cts), 5(cts) and «:(cts) are defined in eq. (5.18), and the prime denotes 
differentiation with respect to their argument. 

In summary, the quadratic approximation to the operator xe Lo(ds, M, N) eq. (5.12) 
is equal to that which is found when cts is an ordinary commuting function, plus the two 
commutator corrections Xi° e<l- (6.19) and xf° eq. (6.21), up to 0{{asPo)‘^) corrections. 
Note that because of the cubic term, the minimum is shifted away from M = 4. For 
instance, the minimum of the symmetrized collinear contributions is shifted to Mq given 
by 


Mo = 


1 

2 


k(as,iV) + Kf^°{aa,N)' 


(6.23) 
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Re-expanding the symmetrized collinear contribntion abont Mq one gets 


M, N) + (M, N) = C(d„ N) + (d„ iV) 


rdo I 


sS{as,Ny 


^ iV) + K^°{as, N) 
+ l{R{as, N) + N)) (M - Mof + 0( (M - Mq)" ) 


(6.24) 


Becanse operator ordering is immaterial in snbleading contribntions, the qnadratic 
approximation at NLO can be obtained by simply adding the leading-order commntator 
corrections Xi° eq. (6.19) and Xs° (6-21) to the NLO kernel x^i^s, M, N) eq. (6.1). 
The qnadratic kernel can then be determined in two eqnivalent ways. One possibility is 
to expand Xai^s, M, N) abont M = |, and determine the shifted minimnm dne to /3o 
corrections and the expansion abont it nsing eq. (6.24), with c and R replaced by the 
intercept and cnrvatnre of the qnadratic expansion of Xai^si N) abont M = i, and 
§00 f^0o replaced by the corresponding parameters of the qnadratic expansion 
of Xi° + Xs°■ This then gives a resnlt of the form of eq. (6.24), which can be pnt on 
shell by letting N = Xo-(cts, -Tf), in tnrn obtained pntting on shell the kernel Xo-(cts, M, N) 
eq. (6.2). A simpler option for the sake of n um erical applications (which we will adopt 
for the compntation of resnmmed anomalons dimensions and splitting fnnctions) is to jnst 
pnt on shell the kernel 

X.(a„ M, N) + alx? (M, N) + xf° (M, iV), (6.25) 


and then determine nnmerically the location of its minimnm Mq and the qnadratic ap¬ 
proximation abont it. Of conrse, the two procednres are eqnivalent np to snbleading 
corrections. 

All the discnssion so far applies to the nnintegrated distribntion, which evolves with 
a symmetric kernel (provided the symmetric operator ordering on the left-hand side of 
eq. (5.14) is adopted). It is easy to switch to the kernel for the nnintegrated distribntion 
by nsing eq. (5.7) in reverse. Specihcally, given the qnadratic approximation 

xli&s^M) = c^{as,N) + \R„{as,N) {M - (6.26) 


to the kernel obtained pntting on-shell the symmetric-variable kernel Xiy{^s,M,N) 
eq. (6.1), switching to the integrated distribntion generates a farther 0{a‘^(3Q) correction 
to the kernel of the form 


with 


xf" = N, M) + N, M) (M - 


cf‘’(d„JV,M) = (M + f) 


-1 


— «s/^o (T7 + y) 


«f°(d„iV,M)= (M+f) 'K(d„iV),M] 

= d;/3 „(A/ + f)-‘A:AA) 


-1 dca(as,N) 
da<i 


(6.27) 


(6.28) 
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where we have taken into acconnt the fact that eq. (5.7) is expressed in asymmetric variables 
so a farther shift is necessary to bring it to symmetric variables as in eq. (5.12). Becanse 
the coefficients and are M-dependent, they mnst be expanded abont M = |. 
Using this expansion in eq. (6.27) and keeping terms np to second order, leads to a set of 
contribntions to the intercept, linear term and cnrvatnre which mnst be added to those of 
eq. ( 6 . 22 ), and treated in the same way in order to obtain the qnadratic approximation to 
the kernel for the integrated distribntion. 

Once the coefficients of the qnadratic expansion abont Mq have been determined, 
they can be nsed to determine the Bateman anomalons dimension eq. (4.23), as well as 
its asymptotic expansion in the /lo —^ 0 limit eqs (4.26) and (4.27). We thns obtain 
the fnll NLO resnmmed anomalons dimension, as discnssed in ref. [7,8] and snmmarized 
in sect. 2 , by determining the pertnrbative rnnning-conpling dnal anomalons dimension 
'Ja NLO, and then combining it with the Bateman resnmmation determined in the qnadratic 
approximation, and snbtracting the donble connting. 

The pertnrbative rnnning conpling dnal of xt,{c(s, M, N) eq. (6.1) is constrncted as 
follows. We start from the observation that [7] the dnal of the NLO donble-leading y 
coincides with the NLO donble-leading 7 , np to snbleading terms. However, recall that in 
order to determine the anomalons dimension we mnst hrst switch back to the integrated 
glnon distribntion. Upon this transformation, Xs changes according to eq. (5.7), while for 
Xo we have 

asXomOiC, M) = asXomMGiC, M) 

= M + 0(d5)]) (?(?, M). 


Now, we note that the dnal anomalons dimension to the kernel Xs(dsM“^)+dsXss(dsM“^) 
which we get after switching back to the integrated level is jnst, by constrnction, the NLO 
anomalons dimension as given by eq. (2.3) np to order On the other hand, the dnal 
of the 0{al) kernel eq. ( 6 . 10 ) is given by eq. ( 2 . 20 ), snpplemented by the 0{al) term 
of eq. (6.29). Note that the resnlt holds with the operator ordering of eq. (6.10). The 
0{as) contribntion to Xa can be bronght to this form by nsing eq. ( 6 . 11 ), and snbtracting 
the donble pole from the /?o term of eq. ( 6 . 11 ), which is a donble-connting contribntion 
between xi and Xs((l ~ M)~^as)- In practice, this means that the /Iq term of eq. (6.11) 
shonld be replaced by Xi° eq. (6.19). 

Pntting everything together, the pertnrbative rnnning conpling dnal to XaNLoiM) 
differs from the hxed-conpling dnal in three respects: hrst, the 0 (ds) contribntion to 
Xe ninst be reordered so that dg is to the left, and second it mnst be snpplemented by 
the contribntion eq. (6.29) which switches to the integrated distribntion; hnally, the hxed 
conpling dnality mnst be corrected by the 0{I3q) term of eq. (2.20). We thns get 


7E''7VLo(«s(t),iV) ^'y^NLo{as{t),N) - Poas{t) 


Xohs{as/N))xohs{as/N)) 

2[x'oMas/NW 


(6.30) 

where the snbtraction in the term in sqnare brackets snbtracts the donble-connting between 
this term and 70 , i.e. it is part of the standard snbtraction cq in eq. (22) of ref. [5]. The 
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anomalous dimension 7 e AfLo (-^7 cts(^)) in eq. (6.30) is the naive (fixed coupling) dual 
eq. (2.9) of Xt.nlo obtained by putting on shell the kernel in asymmetric variables, 


X^nlo{M - + ’ ( 6 - 31 ) 

where Xt.nlo{M, N) is given by eq. (6.1) as discussed in this section, Xi°(-^) is given by 
eq. (6.19), and the term in brackets is the transformation to the kernel for the integrated 
distribution, with the double counting between it and Xs{M~^as) subtracted. 

We finally obtain the full resummed anomalous dimension by combining this result 
with the Bateman running coupling resummation: 

7FiVLo(«s(f), N) = N) - 7 f (a,(t), N) 

JD t:> v ^. 32 J 

- ^ssi(^s{t), N) - Xss,oi(^s{t), N) + 7match(as(t), N) + 7mom(as(t), N). 

In this equation, the Bateman anomalous dimension is given by eq. (4.17), computed with 
the curvature c and intercept k determined from the quadratic expansion of the on-shell 
form of the kernel eq. (6.25) discussed above. However, in oder to ensure complete cancella¬ 
tion of the singularities, the asymptotic double counting terms 7 ^(as(t), N), 7 ^(as(t), N) 
and N) should be determined with values of the parameters c and k which dif¬ 

fer from those used in 7 ® itself by sub-subleading terms. Specifically, 7 ^(as(t), N) should 
be computed with the values of c and k which characterize xt.nlo eq. (6.30), namely, 
the intercept and curvature of the quadratic expansion of the kernel x{1t,nlo) related to 
Xt,nlo by fixed-coupling duality eq. (2.9). These differ from the values of c^° of of 
the quadratic expansion because of the contribution to the latter from eq. (6.21). In 
practice, we eliminate this subsubleading singularity mismatch through the term 


7match (cts (f) 7-^) 






o/3o 




W + 1 


+ 


'N + 1 


+ 


1 + c 


1 -F C^o 


(6.33) 


:K 


^2k{N + 1 ) 


the first two terms shift the square root cut, and the remaining ones ensure that at large 
N the resummed anomalous dimension eq.(6.32) reproduces the NLO GLAP result. Fur¬ 
thermore, 7 ^(as(t), N) is computed with the intercept and curvature of the leading order 
BFKL kernel xo eq. (3.8), which ensures cancellation of the pole at M = ^ of the 0{(3q) 
term in eq. (6.30) which motivates the running coupling resummation, as discussed in 
refs. [7,8] and summarized in sect. 2. The values used in 7 ^ differ from these due to the 
inclusion of all the subleading corrections discussed above, but this is immaterial because 
xfs already subleading. Finally, the term 7 mom serves the purpose of improving the 
matching to large N, by removing terms which survive at the momentum-conservation 
point N = 1. These terms are of course sub-sub leading, so they could be in principle 
omitted. In practice 7 mom is constructed as in eq. (5.21), with dm readjusted to ensure 
that XENLoi^) = 0 - 
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Figure 7: Our most accurate prediction, denoted ‘NLO rc res’, eq. (6.32), is compared to 
the LO, NLO and NNLO GLAP perturbative results and to the DL LO approximation, 
eq. (5.1). 


7. Results for the Anomalous Dimension and the Splitting Fnnction 

We are now finally in a position to collect and discuss our results. Our goal was to 
derive for the singlet anomalous dimension and for the associated splitting function ex¬ 
pressions that at large N or large x reduce to the familiar perturbative GLAP results, and 
at small N or small x are improved by including in a suitably resummed form the informa¬ 
tion from the BFKL kernels with a procedure that takes into full account the constraints 
from momentum conservation and symmetrisation and includes running coupling effects. 
As mentioned in the introduction, all results presented in this section are computed with 
n/ = 0 in order to minimize the impact of issues of scheme dependence. 

The results for the anomalous dimension are presented in figs. 7-8. In fig. 7 our most 
accurate prediction, denoted ‘NLO rc res’, is compared with the LO, NLO and NNLO 
GLAP perturbative results and with the DL LO approximation. The curve labelled ‘NLO 
rc res’ shows the anomalous dimension, which includes the information from the NLO 
BFKL kernel and perturbative anomalous dimension, to which it reduces at large N. It 
corresponds to the result given in eq. (6.32). The DL LO curve, obtained from 'Jdllo 
given in eq. (5.1), corresponds to a naive application of the BFKL kernel xo to evaluate 
the small N behaviour. The step corresponds to the cut starting at the value of N which 
corresponds by duality to the minimum of xo- It is quite evident that our final improved 
anomalous dimension is remarkably close to the GLAP NLO result down to very small 
values of A, where finally we can see a departure due to the position of the Bateman pole 
which is at a small positive value of N and not at A = 0 as is the case for GLAP. The 
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Figure 8: The NLO resummed prediction eq. (6.32), denoted ‘NLO rc res’ is compared 
to the LO and NLO GLAP perturbative results and to the leading approximation, 
eq. (5.20), denoted by ‘LO rc res’. 


instability of the NNLO GLAP result is due to an unresummed ^ pole, and it is removed 
by the resummation. 

Figure 8 shows a magnihcation of the small N region, and also compares the NLO 
resummed result to its LO counterpart eq. (5.20). The stability of the resummed perturba¬ 
tive expansion is manifest in the small difference between the LO and NLO results, which 
however are in perfect agreement with the respective GLAP curves for large N ^ 1. At 
medium-small N the NLO resummed curve is in fact rather closer to the unresummed one 
than the LO resummed. At very small N the LO and NLO resummed curves are almost 
indistinguishable because of the closeness of the respective Bateman poles, displayed in 
Fig. 6, and rise somewhat more steeply than the unresummed result. 

We now consider the splitting function obtained by inverse Mellin transformation from 
these anomalous dimension. In hg. 9 the resummed splitting function curves (labelled 
‘LO rc res’ and ‘NLO rc res’) which are our main predictions at the leading and next-to- 
leading level, corresponding to the anomalous dimensions with the same labels in hg. 8, are 
compared with the LO, NLO and NNLO GLAP perturbative splitting functions (note that 
what is actually plotted is xP). The hrst important feature is the marked deviation of the 
curve NNLO GLAP from the approximate 1/x behaviour at small x which is typical of the 
LO and NLO GLAP singlet splitting function. This pronounced deviation makes the need 
for resummation even more compelling, in the sense that it clearly shows that the ordinary 
perturbative expansion is unstable at small x. As the bulk of HERA data is concentrated 
at X > 10“^ — 10“^ the departure from the NLO splitting functions used in most hts to 
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Figure 9: The LO and NLO resummed splitting function obtained by Mellin transforma¬ 
tion of the anomalous dimensions eq. (5.20) and eq. (6.32) displayed in fig. 8 (labelled ‘rc 
res’) are compared to the LO, NLO and NNLO GLAP perturbative splitting functions 
(note that the plot shows xP). 


the data is not too large, but, for the perturbative results, the difference NNLO-NLO is 
much larger than NLO-LO in the whole range of small x. The resummed NLO prediction, 
which constitutes our main result, closely follows the NLO GLAP curve down to x values 
as small as x ~ 10“^, thus explaining the success of the NLO perturbative QCD fits of 
the HERA data. In comparison the LO resummed prediction shows more pronounced 
deviations from the perturbative splitting function. It is only at very small x, x < 10“®, 
that there starts to be a noticeable difference between the perturbative and the resummed 
splitting functions. This late take-over of the BFKL regime, softened by the resummation 
procedure, is due to the small residue of the Bateman pole. 

In fig. 10 we compare our main results for the splitting function, i.e. the LO and NLO 
resummed curves of fig. 9, to the result obtained in our previous work, ref. [8], denoted by 
‘LO lin rc res’ and with the result (labeled ‘NLO CCSS’) which was obtained by the authors 
of ref. [11] by a different technique based on the same physical ingredients. The prediction 
of ref. [8] is a simple analytical result that can be described as the best approximation that 
can be derived if only the 1-loop kernels 70 and xo are known. Note that, since xo is already 
symmetric for M —> (1 — M), there is no need of symmetrisation, while the result shown 
includes collinear resummation, momentum conservation and running coupling effects. The 
rise at small x is steeper and takes off at larger x than for the more refined results obtained 
in this paper and in ref. [11], because the symmetrised resummed DL kernels displayed in 
figs. 1 and 3 have a smaller value at the minimum and are much flatter than xo- The main 
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Figure 10: The LO and NLO resummed splitting function of fig. 9 are compared to the 
result obtained in our previous work, ref. [8], denoted by ‘LO lin rc res’ and to the result 
of ref. [11], labeled ‘NLO CCSS’ 


advantage of the new resnlt, however, is its perturbative stability, i.e. the closeness of LO 
and NLO results. 

The approach of ref. [11], although based on similar physical principles, differs in 
many respects from ours, most notably in the fact that the whole approach is based on 
the BFKL equation, and attempts to derive the full gluon Green function, instead of 
focussing on the determination of anomalous dimensions, as in our approach, which treats 
the X BFKL kernel and 7 anomalous dimension on the same footing. It is therefore to 
some extent unexpected that our NLO resummed result at small x is very close to the 
NLO CCSS splitting function. In fact, the very small difference between our results and 
those of CCSS is surely accidental, given the uncertainty involved in various subleading 
resummation terms. Note that even though the approach of ref. [11] does not yet allow 
for the inclusion of NLO CLAP terms when n/ 7^ 0, the CCSS curve of hg. 10, taken from 
refs. [12] and computed with n/ = 0, does include a NLO GLAP contribution, which when 
nf = 0 vanishes as x ^ 0. 

A feeling for the uncertainty of our results can be obtained by determining their de¬ 
pendence on the choice scale which in turn may be estimated from the dependence 
of the position of the rightmost singularity Nb of the resummed anomalous dimensions 
eq. (5.20),(6.32) on the value of ag, which was displayed in hg. 6: the value of Nb de¬ 
termines the asymptotic small x behaviour of the splitting function xP ~ We 

X —>^0 

can also use this dependence to estimate the uncertainty on our results due to subleading 
higher order terms. To this purpose, we study the dependence of the value of Nb{c(s) on 
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Figure 11: Dependence of the rightmost pole Nb of the leading and next-to-leading 
order resnmmed anomalons dimensions eq. (5.20),(6.32) on as{kQ^) for various values 
of the renormalization scale k. The solid cnrves correspond to fc = 1 (lower: LO, 
higher: NLO), while the dotted (LO) and dashed (NLO) pair of curves give the range 
of variation as ^ < fc < 4. 


the choice of renormalization scale. At leading order, we plot with Nb the 

rightmost pole of the leading order eq. (5.20), while at next-to-leading order we plot 

NB{k,Q^) = NsiasikQ^)) + Polnkal{Q^)j^NB{as{Q^)), (7.1) 

dcXs 

with Nb the rightmost pole of the next-to-leading order (6.32). In hg. 11 we 

compare the valnes of the leading order and next-to-leading order as a fnnction of 
as{kQ‘^) for /c = 1, /c = 4 and k = j, while in hg. 12 we show the dependence of these 
qnantities on k when is chosen so that as{Q^) = 0.2. We see that the scale dependence 
of the leading-order resnlt is abont as large as the difference between leading and next- 
to-leading order resnlts, whereas the scale dependence of the next-leading-order resnlt is 
smaller by abont a factor fonr. This is somewhat larger than the difference between onr 
NLO result and that of ref. [11], as shown in hg. 10. Also, it appears from hg. 12 that 
whereas at leading order is monotonic as a function of the scale parameter k, at next- 
to-leading order the dependence oi Nb on k is stationary around /c ~ 1, conhrming that 
indeed is the appropriate scale for this process. We conclude from the signihcantly 
reduced dependence on renormalization scale at NLO that we have indeed succeeded in 
constructing a stable resummed perturbative expansion in the small x region. If we were 
to go to NNLO in this expansion, we would expect the scale dependence to be reduced 
still further. 
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Figure 12: Dependence of the rightmost pole Nb of the leading and next-to-leading 
order resummed anomalous dimensions eq. (5.20),(6.32) on the renormalization scale k 
when as{Q^) = 0.2 (lower at large k: LO, higher: NLO). 


8. Conclusion 

In this paper we have presented in detail a complete procednre to constrnct an im¬ 
proved singlet anomalous dimension or splitting function that reduces at large x to the 
NLO perturbative approximation but also contains all resummed small x improvements 
from NLO BFKL kernels. The problem of explaining the apparent smallness of the small 
X corrections in the HERA data region in spite of the wild behaviour of the BFKL pertur¬ 
bative expansion is solved. The main physical reason for the moderate impact of the small 
X corrections is the duality relation that implies that in the large domain where both 
and 1/x are large the GLAP and the BFKL expansions are both valid and know about 
each other. Thus the corresponding leading twist terms must match and indeed a better 
convergence is obtained by the ‘double leading (DL)’ expansion where at each level the 
BFKL (GLAP) perturbative terms are used to resum the corrections to the GLAP (BFKL) 
kernel in a symmetric way. The constraints from momentum conservation and from the 
underlying symmetry of the gluon-gluon BFKL scattering kernel under interchange of the 
two gluons are used to construct a stable and physically motivated expansion at each level. 
An important effect is also obtained from the running coupling corrections which are large 
and considerably soften the small x asymptotic behaviour. The technical implementation 
of the running coupling corrections in our approach is realized via a quadratic approxi¬ 
mation of the BFKL kernel in the DL approximation near its central minimum followed 
by the analytic solution of the corresponding differential equation in terms of a Bateman 
function. The systematic use of duality and the Bateman method are the main qualitative 
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differences with respect to the approach of ref. [11] which uses a different path to reach 
similar conclusions from the same basic physical principles. 

The main results of our work for the singlet anomalous dimension and the correspond¬ 
ing splitting function are displayed in hgs. 7-11 where the comparison with the results of 
ref. [11] is also shown. In the region of the bulk of HERA data for x > 10“^ the corrected 
splitting function at NLO level shows a moderate depletion with respect to the pertur¬ 
bative NLO approximation. At still smaller values of x the asymptotic regime eventually 
takes over and a powerlike increase of the splitting function is observed, though at a much 
lower rate than that naively expectated from the value of the Lipatov exponent. 

We recall once again that to make direct connection to the structure functions we 
need to include the coefficient functions. This in principle is relatively easy, because the 
coefficient functions are available, even though properly accounting for scheme choices, 
diagonalization of the anomalous dimension matrix, and matching to large x coefficient 
functions is a laborious task, as discussed in ref. [6]. The effects of the coefficient func¬ 
tions are subdominant in the factorization scheme adopted here, but they deserve a full 
investigation. 

In conclusion we think that by now the problem of the small x behaviour of the singlet 
splitting functions has been well understood in its physical principles and that the technical 
tools have been satisfactorily developed. The next step is to make a direct comparison of 
the resulting theoretical predictions with the data. 
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Scottish Universities Physics Alliance. 


Appendix A 

To evaluate xi(M) we use the result 


Jo ^ ~r X JX ^ 

TT^ TT^ 

GsIuttM 6M(1 —M) 


where 


and 


OO 


r=l 




(-)' 


—^ {m + k + 

k=l m=l ^ A' 


(A.2) 


(A.3) 
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so the sum over r has inhnite radius of convergence. If we further note that 


sin ttM 


= + f) - ) + (M ^ 1 - M), 


(A.4) 


cos ttM 
sin^ ttM 




(A.5) 


then we may then use eq. (14) of ref. [16] to write xi(M) = -?^5{M) — as 


XI (M) 


-i/^o^(^Xo(M)2 - - V^'(l - M)) 

+ A [(f - X - ^>(’^(1) - >^0'/))+3C(3)+ r(M) 

+ 4(ir^(V>(^ + f-) -’/'(-f)) + ^ + “ J»>.(V)) 


4(1 -2M) 


3 + (l + ^) 


(2 + 3M(l-M)) 
(3-2M)(l + 2M) 


)(V^'(i + f)-V^'(f)) 


+ (M ^ 1 - M) . 

(A.6) 

The advantage of using this expression is that all terms in square brackets are regular for 
ReM > 4 (note that the poles at M = 4^ M = | and M = — 4 in the last term are 
all removable: they cancel against similar poles in the M —1 — M piece). Thus we may 
extend them off-shell following the same procedure as at LO. The hrst term may be treated 
by the simple observation that up to subleading terms we may substitute agXoiM) = N, 
which is trivially extended off-shell as N = asXo{M, N). The hnal result for xi(Af, N) is 
then (in asymmetric variables for simplicity) 


Xi(M - f, iV) = -i/3o^(^Xo(M - f, Nf - i;\M) - - M + N)) 

+ S[(f-^- + N)) - V^(M)} + 3C(3) + r{M) 

+ 4{7r^(V'(^ + f)- + ^ + lk(^(l) - ^(1 + ^)) - 

+ + f) - nf)+^\i) - v^'d)) 

3M(1 - M))| + T") ~ ~ '*A^(|)) 

+ 2(T^(d(i + f) - d(f) + d(-i) - d(i)) 

+ 2i3-2M)('^'(l + ¥) 

+ ^(1 + - V^'(l) + 3d(|) - 3V^'(4))iV 

+ (M ^ 1 - M + iV) . 

(A.7) 

where in the term in curly brackets we have taken care to ensure that we do not spoil the 
cancellation of the poles at M = ±4, M = | by adding in subleading (ie 0{N)) terms: 
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the growth of these extra terms at large N must then be subtracted off again, and this is 
done in the last line. The result in symmetric variables may be found by the simple change 
of variables M —M + ^. 

The off-shell extension of eq. (6.14) is then straightforward: 

Xi(M, N) = xi(M, N) - ixo(M, N)^ [2^P\1) - + f) - V^'(l - M + f)] . (A.8) 

The collinear subtraction eqn(6.15) is then 


Xi(M,iV)=xi(M,iV)-^i 


M+ 


W 


+ 


1-M+- 


W ~ 92 


+ 


(«+f) 


(A.9) 


with coefficients gi as before, eq. (6.16), which render in particular Xi{M — ^,iV) hnite 
in the collinear limit M —0. To improve the matching to GLAP at large N, we may 
further add to Xi(-^ — A^) a manifestly subleading contribution — N) -t-Xi(0, 0) 

to remove spurious constant terms at large N, and thus spurious 0(1/A^) terms in the 
anomalous dimension. 
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